Linear response theory for quantum open systems 
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Basing on the theory of Feynman's influence functional and its hierarchical equations of motion, 
we develop a linear response theory for quantum open systems. Our theory provides an effective 
way to calculate dynamical observables of a quantum open system at its steady-state, which can be 
applied to various flelds of non-equilibrium condensed matter physics. 



Introduction: The linear response theory (LRT) has 
been widely used in condensed matter physics since it was 
derived by Kubo in 1957 [1]. For example, one can calcu- 
late the conductivity via the current-current correlation 
function directly from the ground state without apply- 
ing any realistic bias voltage. However, Kubo's theory is 
only valid for an equilibrium closed system (not limited to 
quantum one). In recent years, quantum theory for open 
systems have attract more and more research interests in 
molecular electronics, nanophysics and biophysics, etc. 
One fundamental issue is whether we can directly calcu- 
late the dynamical observables (e.g. spectra function) of 
a static open system (with time-translation-symmetry) 
instead of evolving it under some time-dependent field? 
In this letter, we address this issue by developing a lin- 
ear response theory for quantum open systems (in anal- 
ogy with Kubo's). Our derivation is based on the theory 
of Feynman's influence functional |2] and its hierarchical 
equations of motion (HEOM) [3]. 

There exist two correlated questions in developing the 
LRT for quantum open systems, 1) how to get the re- 
duced density operator (RDO) of a static open system; 
and 2) how to calculate the response to another probe 
field? Since question 1 has been well resolved recently 
via the HEOM of RDO (and its auxiliary ones) derived 
by differentiating Feynman's influence functional we 
thus focus on question 2 in the present work. 

Before the formal derivation, let us introduce an im- 
portant character of influence functional: if a system S 
simultaneously interacts with two different environments 
(A and B) and no direct coupling between them exists at 
initial conditions, then the total influence functional F of 
the system is: F = Fa ■ Fb 0- It also means the influ- 
ence phase satisfying: $(i) = $^(0 + <i>s(i) [2]. When 
one differentiates the influence functional to derive the 
HEOM fs'l, those two different phases will linearly enter 
into HEOM tie by tie with their simple relationship be- 
ing maintained [see Eq.([5])]. That character provides a 
practical solution for question 2 mentioned above. 

HEOM of propagator: As the answer of question 1 
(and for later purposes), we first outline the main re- 
sults of RefH in the language of Liouville-space propa- 



gator. Let us suppose our quantum open system com- 
posed of the (reduced) system with particle operator aj^ 
(tj = -\- / — corresponding to the creation/annihilation op- 
erator ) and the environment with bath correlation func- 
tion C^^{t) (For fermion bath, cr = +/— corresponding 
to particle transferring in/out of the system), where /i {v) 
represents the orbital (site, energy or spin, etc.) index of 
the system. We then chose {ip} as an arbitrary basis set 
(defining a certain path in path integral representation) 
in the system subspace, and = {tp,ip'} for short of the 
two paths in the theory of influence functional [2J. Now 
we can write the HEOM of the influence functional in 
Liouville space as [j = {fia}, j = {//tr}] 
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where the sum runs over a\\ j ^ jk] k — 1, ...,n. F^^"^ 
is the superoperator form of the ordinary influence func- 
tional F, and the auxiliary ones. The expressions 
of auxiliary influence functionals appearing in the rhs of 
Eq.H]) are Q 
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(2c) 



The auxiliary RDOs, , can hen be deflned in terms 

(3) 



of the auxiliary influence functionals 

pf\t)^ul''\tM)p{to). 
where Uj"'\t,to) is the time-evolution superoperator 

'^o[*o] 



with S[ip] being the classical action functional of the sys- 
tem. 

From Eq.dT]), one can derive the HEOM of RDO as 
shown in Ref.'Si. Here, we give the HEOM of the reduced 
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Liouville-space propagator, defined as ^/j (t — to) = 
Z^j"^(i,to) for the time-translation- invariant system 

g{t) = -[iCs + ic,f{t)]g{t) -ij:^ A-^g^^\t); (5a) 
x(-i)"-^-4eir'^W -*E/A^e;;+'^(t). (5b) 

with the initial condition being ^?j"''(i = 0) = 6no- In 
EqlSl £s is the Liouville operator of the system, £s--- = 
[Hs, ...], while Csf(t) is that of arbitrary external time- 
dependent field (e.g. the probe field in the LRT). 

HEOM space : We now define the HEOM linear space 
that can be seen as an extension of Liouville space. In 
the former the basic element is no longer a operators as 
in the latter, but is {N -I- 1) dimensional super-vector 
constituted by the operator and its auxiliary ones. Let 
us take the RDO as an example, which is extended as 
the HEOM-space RDO, i.e. 

p{t)^{p{t), pf{t)[n^l,2...N]). (6) 

According to Eqs.([3]) and (0), the time-evolution of p{t) 
is determined by the HEOM-space reduced propogator 

Q{tM) 

p{t) = g{t,to)p{to). (7) 
In HEOM space, Eq.([5]) can be shorted as 

dg{t,to)/dt = -Ait)git,to), (8) 

where A{t) is the main time-evolution superoperator act- 
ing on super-vectors in HEOM space, which determines 
the equation of motion (EOM) of the propagator as well 
as the density operator 



N 

{{A\B)) ^ {{A\B)) + ^ ((Af^li?]"^)) , (11) 

ri=l 

where {{A\B)) = Ti-[A+B] being the standard defi- 
nition of inner product in Liouville space. Similarly, 

(«>|i,<-)).T.{[4f B<"'}. 

LRT in HEOM space: We are now on the position 
to develop the linear response theory in HEOM space. 
Suppose a variable A of the system S is perturbed by 
a weak probe field epr{t), which linearly couples to the 
system via its another variable B 

Hpr{t) = ~Bepr{t). (12) 

For simplicity, we suppose both A and B are Hermitian 
operators (extending to non-Hermitian cases is straight- 
forward) . The change of the expected value of A caused 
by the disturbance of epr(t) is 

6A{t) = Tr [ASp{t)] = {{AmSpm , (13) 

where A{0) and 6p{t) respectively denote the expansion 
of A and Sp{t) in HEOM space, i.e. 

A(0) = {A, 0}; (14a) 
5p{t) = {(5p(t), (5pj"^(n = 1,2...7V)}. (14b) 

We then apply the first order perturbation in the EOM 
of RDO, 

p{t) = p,{t) + Sp{t) = g{t,T)p{T). (15) 

By separating the superoperator A{t) as A{t) = Ag + 
Apr{t), and inserting the Dyson's equation in HEOM 
space [tI into EqfTSl 



Pit) = -Ait)p{t). (9) 

Although the concrete form of A{t) can not be written 
out separately from Eq.(l8]), one can easily work out its 
effect from Eq.® 

Generalizing above definitions, any operator in Hilbert 
space now can be expanded to a super-vector in HEOM 
space, which obeys similar EOM as Eq.®, 

A{t)^{A, A^\n^l,2...N)}. (10) 

The inner product of two vectors A and B in HEOM 
space is defined as 



g{t, r) = gsit -r)- J dT'g{t,T')Apr{T')g,{T' - r), 

(16) 

we arrive at 

Sp{t)^-f dTgit,T)Apr{T)p,{T). (17) 

In above derivation, we have used ps{t) = ^s(i — t)Pj,(t) 
and initial conditions Sp{T = 0) = 0; p{t = 0) = 
P.(t-O). 

To proceed, we define the time-independent HEOM- 
space superoperator S as 
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B = lApr{t)/epr{t), (18) 

whose action can be determined from Apr{t)A — 
i[Hpr{t), A] [see Eq.®] as 



Bp=[B,p]. (19) 
Inserting Eg. p^ into pT|). we finally get 



Spectra function in HEOM space: In principle, one can 
calculate the dynamical correlation of any two system- 
operators {A and B) via the HEOM-space LRT. One 
typical example is the spectra function that plays im- 
portant role in many body physics. In what follows, 
we will demonstrate how to obtain the (reduced) spec- 
tra function of a quantum open system. The spectra 
function Jab{^) directly relates to the imaginary part of 
the retarded single-particle Green's function G^g^t) in 
the form of 



SAit)^t dT((A(0)|g(t,r)B|p(T)))ep,(r), (20) 

from which we are ready to define the response function 
in HEOM space as 



XAB{t,T)=t{(A{0)\git,T)B\p{T) 



(21) 



Obviously, the definition of Eq. (PT|) is quite gen- 
eral. For a static quantum open system satisfying time- 
translation-symmetry, p{t) — p^^{T) = {peq, Peq^[n, = 

1, 2...iV]} and g{t, r) -J> Q^it ~ r), which leads to follow- 
ing definition of the response function 



XAB{t) 



A{0)\gs{t)B\p,JT) 



(22) 



The correlation function in HEOM space can be simi- 
larly defined as |8|] 



CAB{t)^((AmgS)B\p,^{T) 



(23) 



where Bp^g{T) = {Bp,g, Bp%\n = 1,2...A^)}. 

For practical calculation, we can rewrite Eq. (|22p as 



G^ABit) ^ -ie{t){{A{t), B})- (27a) 
Jab{lo) = --Iu,[G^ab{^)]- (27b) 

TT 

Please be noted that the Green's function in our theory 
is defined for two arbitrary operators, which can reduce 
to the one in textbooks by setting A — a and B = . 

For fermion, G'^^(t) can not be obtained directly from 
XABit) due to the anti-communication relation in the 
former but communication one in the latter. Fortunately, 
we can get G^^(i) from correlation function CABit), i.e. 



G'ABit) = -^0{t){{A{t), B}) 

= -i0it) \CAB{t) + CBA{-t) 



(28) 



To proceed, we introduce the general spectra func- 
tion Cab{^) = \ dt e^'^^CABit), which satisfies the 
detailed-balance-relation Cba{~'^) ~ e '^"C'ab('^)- Af- 
ter some algebra, we obtain 



JAB{io)^-{l + e-^^)CAB{oj), 



(29) 



which is obviously the fluctuation-dissipation theorem in 
HEOM space. 

Since in Ea. (|23p only ^^(0 is the function of time, we 
have 



XABit) = ((A(O)k(i))) = Tr [A+a{t)] , (24) 



where 



Cab{(^) 



AmOs{to)B\p^g{T))), (30) 

Thus 



where 5s (w) is the Fourier transform of Eq 
JABi^) can be calculated from 



(Tit) = g.(t)*T(0); (25a) 
aiO) = *Bp,,(T) 

= {t[B, p,,], t[B, p(^)](n = l,2...7V)}.(25b) 

Similarly, if we set cr(0) = Bpg^(T), Eq.lHSD can be 
rewritten as an easier handling form 

CABit) = ((A(O)l^(i))) = Tr [A+ait)] . (26) 



JABi^) = ^{l + e-^-) 



AiO)\g4uj)B\p,JT) 



(31) 



More concretely, if choosing cr(0) = Bp^^iT) 
{Bpeq, Bp''e^\n — 1,2...A^)}, then we have cr(w) 
as(w)o-(0) and 



JABiu^) = ^(l-He-^-)((A(0)|<T(c.))) 
= ^(l + e-^")Tr[A+a(..)]. 



(32) 
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Eg. (13211 is the main formula to calculate HEOM-space 
spectra function in our theory. Before that, one must 
solve equation Agp^^iT) = to obtain p^^{T). 

Summary: In summary, we have developed a linear re- 
sponse theory for quantum open systems on the basis of 
the theory of Feynman's influence functional and its hi- 
erarchical equations of motion. From our theory, one can 
directly calculate the dynamical observables (e.g. spec- 
tra function) of a static open system instead of evolving 
it under time-dependent field. It can be applied to non- 
equilibrium many-body physics, nanophysics, etc. 
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